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The energy dependence of particle production sources in relativistic heavy-ion collisions is in¬ 
vestigated from RHIC to LHC energies. Whereas charged-hadron production in the fragmentation 
sources follows a ln(sA^iv/so) law, particle production in the mid-rapidity gluon-gluon source ex¬ 
hibits a much stronger dependence oc \n^(sNN/so), and becomes dominant between RHIC and 
LHC energies. The production of particles with pseudorapidities beyond the beam rapidity is also 
discussed. 
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I. INTRODUCTION 

The investigation of charged-hadron production in rel¬ 
ativistic heavy-ion collisions has generated a vast amount 
of energy- and centrality-dependent data at energies 
reached at both, the Relativistic Heavy-Ion Collider 
RHIC [H , and the Large Hadron Collider LHC [H . It has 
been shown Q within the framework of a nonequilibrium- 
statistical relativistic diffusion model (RDM) Hi that 
the energy-dependent multiplicity of produced charged 
hadrons is well understood quantitatively based on a 
mid-rapidity low-x gluonic source and the two fragmen¬ 
tation sources. This applies not only to An An collisions 
at RHIC [il and PbPb at LHC [ 3 ] , but also to asymmet¬ 
ric systems such as dAu at RHIC [H and pPb at LHC 

i- 

The relativistic diffusion model is in scope and charac¬ 
ter located between the (equilibrium) statistical model 
for multiple hadron production that was proposed by 
Fermi [9[ and Hagedorn [T^ , and much more detailed nu¬ 
merical models that aim at a microscopic description of 
the collision, such as the Color Glass Condensate (CGC, 
see 0) for the initial state, hydrodynamics for the main 
part of the time evolution (e^. |12l-[l^). and codes like 
URQMD for the final state p^ . 

The statistical hadronization (or thermal) model has 
been further developed and compared to a large amount 
of data by many authors such as Braun-Munzinger et 
al. or Becattini et al. [iM3, and it has consistently - 
with only few exceptions - provided good descriptions of 
particle production yields, in particular, at mid rapidity. 
As a consequence of its ambition to account for particle 
production with few parameters (temperature, chemical 
potential, characteristic volume) in an equilibrium set¬ 
ting with collective expansion, the thermal model does, 
however, not describe effects such as the plateau occur¬ 
ring in rapidity distributions dN/dy of produced particles 
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at higher (RHIG and above) energies, the corresponding 
dip in pseudo rapidity dN/dr], and other outstanding fea¬ 
tures such as limiting fragmentation at RHIG and LHG 
energies. 

To account for such non-equilibrium effects and model 
the collision in full detail requires in current scenarios 
to match the GGG initial state smoothly to viscous hy¬ 
drodynamics when the coupling constant becomes too 
strong in the course of the time development for pertur¬ 
bative QGD techniques to be applicable [2^, and finally 
use Gooper-Frye freeze out HH or another code that ac¬ 
counts for the final-state interactions [T6|. 

However, since even the most sophisticated codes that 
purport to describe the full time evolution will contain 
a certain amount of arbitrariness and can not fully re¬ 
place the experiment, it appears indicated to permit phe¬ 
nomenological models such as the RDM that include non¬ 
equilibrium effects to some extent, reproduce substantial 
features of the data and have some predictive power, but 
do not claim to fully account for every detail of the col¬ 
lision and of the ensuing particle production. 

The nonequilibrium-statistical relativistic diffusion 
model is - in its linear approximation (Ef - based on 
an analytically solvable transport equation with three 
sources. It does not only consider particle production 
from a central source as the thermal model does, but 
also from the fragmentation sources. The latter evolve in 
time and eventually tend to merge with the central source 
towards an overall thermal equilibrium distribution, but 
since the interaction time is extremely short at RHIG 
and LHG energies, this equilibrium state is not reached, 
and in particular the rapidity and pseudorapidity distri¬ 
butions show characteristic nonequilibrium features. 

In this work I present an investigation of the energy 
dependence of the charged-hadron production sources 
within the relativistic diffusion model in symmetric sys¬ 
tems, An An at RHIG c.m. energies per nucleon pair of 
19.6, 62.4, 130 and 200 GeV, and PbPb at LHG ener¬ 
gies of 2.76 and 5.52 TeV. The gluon-dominated source, 
in addition to the fragmentation sources related to the 
valence part of the nucleons, had been implemented ear- 
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Her into the RDM Hi- A related model with a gluonic 
source at mid rapidity had also been proposed by Bialas 
and Czyz [22|. 

In 0 it has been found that the fragmentation sources 
for produced charged hadrons - which are clearly visi¬ 
ble in net-proton rapidity distributions where the gluonic 
source cancels out [23| - have the expected logarithmic 
dependence on ^/sNN^ whereas the particle content in 
the mid rapidity gluon-gluon induced source that rises 
strongly with energy is close to a power law. This result 
has since been corroborated through other independent 
investigations of charged-particle and transverse energy 
production [ 23 , [HI such that a renewed and more pre¬ 
cise consideration in particular of the central source is 
indicated. 

The fragmentation sources are responsible for most of 
the yield in the regions close to the beam rapidities. Here 
limiting fragmentation scaling [h is valid not only at 
RHIC, but also at LHC energies Q- This is in contrast to 
earlier predictions of the thermal model [HI which hud a 
violation of extended longitudinal scaling at LHC ener¬ 
gies, providing another indication that equilibrium sta¬ 
tistical concepts are invalid in the fragmentation region. 
Here the yields in pseudorapidity also extend beyond the 
value of the beam rapidity, and in the hnal paragraph of 
this note the origin of this effect is discussed. 



n 

FIG. 1. (Color online) The RDM pseudorapidity distribution 
functions for charged hadrons in central AuAu (RHIC) and 
PbPb (LHC) collisions at c.m. energies of 19.6 GeV, 130 GeV, 
200 GeV, 2.76 TeV and 5.02 TeV shown here are optimized 
in hts with respect to the PHOBOS 0,113 (bottom) and 

ALICE [^ (top) data, with parameters from [^. The upper 
distribution function is an extrapolation to the LHC design 
energy of 5.52 TeV. At the lowest energy, only the fragmen¬ 
tation sources contribute (dash-dotted curves). 


II. HADRON PRODUCTION SOURCES 


For a detailed phenomenological investigation of the 
charged-hadron particle content in the three particle- 
production sources, the nonequilibrium-statistical rela¬ 
tivistic diffusion model is used. The fragmentation 
sources Ri^ 2 {y^t = '^int) with charged-particle content 
(projectile-like), (target-like) and the midra¬ 

pidity low-x gluon-gluon source Rgg{yR = Tint) with 
charged-particle content are added incoherently to 
generate the total pserudorapidity density distribution as 


dNl^{y,t = Tint) 
dy 


N^J’^R,{y,nnt) 


+N^J’^R2{y, Tint) + Rut) 


( 1 ) 


with the rapidity y = 0.5 • ln((F^ p)/{E — p)), and the 
interaction time Tint- The latter corresponds to the to¬ 
tal integration time of the underlying partial differential 
equation, which is a linear partial differential equation of 
the Fokker-Planck type, as described in Q. 

Converting the rapidity distribution dN/dy for pro¬ 
duced charged hadrons to the corresponding pseudora¬ 
pidity distribution dNjdr] {f] = - ln(tan(6>/2)) ) with the 
proper Jacobian transformation dy/dr] and minimizing 
the analytical solutions of the transport equation with 
respect to available pseudorapidity data then yields the 
particle content of the sources as functions of ^JsNN 0. 
The corresponding RDM-parameters for central collisions 
have been published in Tab. 1 of 0. 


Results of this approach are summarized in Fig. [U 
where the charged-hadron pseudorapidity distributions 
are shown from low RHIC energies of 19.6 GeV, via 130 
GeV, 200 GeV, to 2.76 TeV, plus a prediction at 5.52 
TeV. It is noted that the midrapidity source is found to 
be absent at 19.6 GeV and appears only at the higher 
energies, rising in particle content with s/snn- The in¬ 
dividual sources are displayed in Fig. [2] at 200 GeV and 
2.76 TeV, where the effect of the Jacobian transformation 
from rapidity y to pseudorapidity rj is also shown. The 
central gluon-gluon source is seen to become dominant 
as the energy is increased from RHIC to LHC. 

The corresponding particle contents of the sources are 
displayed in Fig. [3l which resembles the analogous hgure 
in 0 , but differs in a decisive detail. The total particle 
content is found to follow a power law, 

= 1-1 • 104(siviv/so)°-'" (2) 

with So = 1 TeV^, whereas the particle content in the 
two fragmentation sources is as expected a logarithmic 
function of the energy 

= 695 . Hsnn/so) (3) 

with So = 100 GeV^. The midrapidity gluon-gluon source 
is approximated by the thin dashed line following a power 
law as was already proposed in 0 

TVf-4.1o3(sxivAof^^ (4) 

with So = 1 TeV^. However, when considering also 
the yield predicted within the relativistic diffusion model 
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FIG. 2. (Color online) The RDM pseudorapidity distribution 
functions for charged hadrons in central 200 GeV AuAu (top 
frame) and 2.76 TeV PbPb collisions are adjusted through 
minimizations to the PHOBOS [l[ (see also [^) and AL¬ 
ICE [ 2 ^ data, see i. The underlying particle production 
sources are shown: dash-dotted curves are the fragmentation 
sources, dashed curves the mid rapidity gluon-gluon sources, 
and dotted curves the central sources without the effect of 
the Jacobian transformation from rapidity to pseudorapidity. 
The particle content in the gluon-gluon source rises strongly 
with increasing c.m. energy, and constitutes the largest source 
at LHC energies. 


(RDM) for the LHC design energy of 5.52 TeV, the power 
law fails to fit the expected yield, whereas a cubic log de¬ 
pendence agrees with the prediction, 

Vf® = 7.5 ■ Inysjvjv/so) (5) 

where 5o =169 GeV^. 

It remains to be seen whether the data actually fol¬ 
low the model prediction. In the upcoming PbPb run 
at the LHC in 2015, the c.m. energy is scheduled to 
be 5.125 TeV, corresponding to 13 TeV pp. The total 
charged-hadron yield predicted by Eq. (2) at this energy 
is = 23,327, with the central source contributing 
^ch ~ 12,811 charged hadrons according to Eq. (5). 
The RQM-value for the total charged-hadron production 
at the lower LHC ener^ of 2.76 TeV is = 17,327 
according to Tab. 1 of Q; the power law Eq. (2) yields 
17,546. The ALICE collaboration meanwhile quotes an 
extrapolated value of 17,146 ± 722 j^ . 

When examining the RDM results for the particle con¬ 
tent of the sources more closely also in the low-energy 
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EIG. 3. (Color online) Number of produced charged hadrons 
as function of the c.m. energy ^'snn from RDM-fits of the 
available data for central heavy-ion collisions at 0.019, 0.062, 
0.13, 0.2 TeV (RHIC, AuAu), 2.76 TeV (LHC, PbPb), plus 
extrapolation to 5.52 TeV. Circles are the total numbers, fol¬ 
lowing a power law oc Triangles are particles from the 

fragmentation sources oc log(sww/so)- Squares are hadrons 
produced from the midrapidity source, with a dependence 
oc \og^{ snn/ so). A power law oc [ 3 ] is also shown (short- 
dashed curve), but fails to fit the extrapolated 5.52 TeV yield. 
The gluon-gluon source (dashed) becomes the main source of 
particle production between RHIC and LHC energies. 



EIG. 4. (Color online) The total charged-hadron produc¬ 
tion in central AuAu and PbPb collision in the energy re¬ 
gion 19.6 GeV to 5.52 TeV is following a power law Ntot oc 
{snn/ so)^'‘^^ (solid line), whereas the particle content in the 
fragmentation sources is Nqg oc In (satat/sq), dash-dotted 
curve. The particle content in the mid-rapidity source obeys 
Ngg OC In^ {snn/so), dashed curve, not too far from a power 
law (short-dashed line) only in the intermediate energy range 
0.1-2.76 TeV. The energy dependence of the mid rapidity 
yield is shown as a dotted line, with PHOBOS data [I|] at 
RHIC energies, and ALICE data [23] at 2.76 TeV. 
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region where RHIC data are available, it turns out that 
the power law Eq. (4) is an acceptable approximation to 
only between about 100 GeV and 2.76 TeV. 

This becomes particularly obvious in Fig.lH where the 
same plot is shown using a double-logarithmic scale, fol¬ 
lowing a suggestion by Trainor [3l|. Here power laws 
appear as straight lines - such as the one for the total 
charged-hadron production, or also for the midrapidity 
yield 


—= 1-15 • 10^{snn/so)^'^^^ ( 6 ) 

with So = 1 TeV^ (dotted line, and data points from 
Phobos [H and ALICE (H). 

The cubic-log dependence of the gluon-gluon source 
(dashed) is seen to fit the points extracted from the 
RDM-analyses Q of PHOBOS and ALICE data rather 
precisely at the available energies, and it agrees with the 
RDM-prediction at the LHC design energy of 5.52 TeV. 

As required by the RDM analysis of the 19.6 GeV 
AuAu data, the gluon-gluon contribution becomes unim¬ 
portant below 20 GeV - whereas a power law would 
still predict a yield of about 100 charged hadrons in this 
energy region. Although a hybrid function with a log- 
dependence at RHIC energies that turns into a power law 
at LHC energies may appear as a reasonable compromise 
nisi, it can not compete with the In^-dependence for 
the central source regarding the precision of the fit to the 
RDM-results. 


III. ENERGY DEPENDENCE OF THE 
MID-RAPIDITY SOURCE 

The origin of the cubic-log dependence of the total 
charged-hadron yield on snn (or -\/snn) in the mid¬ 
rapidity gluon-gluon source can be traced schematically 
neglecting for the moment the precise value of the pro¬ 
portionality factor appearing in Eq. (5). The width of 
the gluon-gluon distribution is expected to scale roughly 
with the beam rapidity, 

o- cx j/beam = In (V^NN/rtip) = 0.5 In {sNN/ml) (7) 

where rup is the proton mass. With respect to the midra¬ 
pidity value, the STAR collaboration observed in 2004- 
2006 that dijet production which generates the hard 
component of the spectrum is at midrapidity propor¬ 
tional to the square of the soft-component density, that 
is associated with \ow-x gluons [32, • 

Since the yield of \ow-x gluons is proportional to the 
logarithm of the c.m. energy, the density at midrapidity 
that arises from dijet production is proportional to In^ s. 
Hence, the integrated yield in the gluon-gluon source can 
be estimated as 


/ yheam ^ AT 

— \ggdr)(x\v?{sNN/so) 

■?/beam ' 


(8) 


in agreement with the above result of the phenomenolog¬ 
ical RDM-analysis. 

On the theoretical side, the gg gg scattering am¬ 
plitude has been evaluated in the presence of a classical 
color field e.g. by Cheung and Chiu Q. They find 
that the classical color field modifies the gg ^ gg elastic 
scattering amplitude, and suppresses it when the longi¬ 
tudinal momentum fraction x of the incident gluon is 
small. The rise of the cross section with energy in the 
central distribution - that is driven by the growth of 
the gluon density at small x - is therefore suppressed by 
the quantum-classical interaction from the dense medium 
Q. The predicted cross section has a In^s asymyptotic 
behavior that satisfies the Froissart bound [35|, and the 
integral over rapidity becomes proportional to In^ s. 

It is interesting to compare the results of the present 
analysis with the rapidity distributions from the hydro- 
dynamic approach of Landau and Belen’kji , and 

applications to particle production by Carruthers and 
Duong-van [sl, |39| , as well as Steinberg [40| . There the 
width (FWHM) T = V8\K2 •cr of a gaussian pseudorapid¬ 
ity distribution for produced charged particles is obtained 
from the variance j39| 

<^Landau = In 7 = In (yiVv/2mp) (9) 


with the Lorentz-factor 7 = l/y^l — /d^),/3 = p/E. It 
turns out that this expression is in reasonable agreement 
0 with data from AGS and SPS where the stopping 
fraction is sizeable. Deviations start to become visible 
at RHIC - where nuclear transparency with well- 
separated fragmentation sources is already obvious - and, 
in particular, at LHC where the measured width of the 
dN /^^-distributions for charged hadrons is substantially 
broader than predicted by Eq. m, as shown by the AL¬ 
ICE collaboration [^. It has therefore been concluded 
’...that Landau hydrodynamics does not explain the ex¬ 
pansion dynamics at LHC energies ’ M- 

Whereas this is certainly true for the overall pseudo¬ 
rapidity distribution of charged particles, Landau’s ap¬ 
proach may still be viable for a proper description of the 
mid-rapidity source which accounts for particles gener¬ 
ated from \ow-x gluons. Indeed for 2.76 TeV PbPb, the 
RDM-analysis yields a width in rapidity y of = 6.24 
0, compared to a Landau result of = 6.36 in 

T] according to Eq. At RHIC energies, the Landau 
result is, however, larger than the RDM result for the 
mid-rapidity source, and the results at the higher LHC 
energies of 5.125 TeV and 5.519 TeV PbPb remain to be 
seen. 


IV. YIELDS BEYOND THE BEAM RAPIDITY 

Already in the investigation of AuAu collisions at 
RHIC energies [l| it has been observed that pseudora¬ 
pidity yields of produced charged particles extend sig¬ 
nificantly beyond the value of the beam rapidity. This 
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is particularly obvious in PHOBOS AuAu results at 130 
GeV where dN/dr] data have been taken beyond ^beam 
ji^ . The RDM solutions for 200 GeV AuAu and 2.76 
TeV PbPb also clearly indicate expected yields beyond 
Vheam at these higher energies, see Fig. [5l 

Obviously it is not excluded that this can partly be 
due to a real physical effect, with a few charged particles 
produced at larger rapidities than that of the beam value. 
However, the bulk of the large charged-particle pseudo¬ 
rapidity density in the region at and beyond the beam 
rapidity - which amounts to more than 100 charged par¬ 
ticles - is likely due to the transformation from rapidity 
to pseudorapidity. 

Reconsider the expressions for rapidity y, longitudinal 
velocity /3||, and pseudorapidity 77 



n 



1 + / 3 || 

1 -/3|| 


exp {2y) - 1 
exp (2y) + 1 


/-|Q\ FIG. 5. (Color online) Produced charged particles in central 
^ AuAu collisions at ^ysNN = 130 and 200 GeV (RHIC/ PHO¬ 

BOS data [I[, bottom, and in 2.76 TeV PbPb Q, top, in com¬ 
parison with the RDM solutions. The values of the beam ra¬ 
pidities are indicated as arrows (ybeam = 4.932, 5.362, 7.987). 
(11) The pseudorapidity yields extend beyond ybeam, which is par¬ 
ticularly evident in case of the 130 GeV PHOBOS data. 


77 = — In (tan( 6 >/ 2 )). ( 12 ) 

The transformation between 77 and y is 


I \ (m/pT)‘^ + cosh^ y + sinh 77 

y = P'^~i= =-• 

W (m/pTY + cosm y — sinh 77 

Here m is the mass of the particle species considered. 
The relative particle abundances in central (0-5%) PbPb 
collisions at 2.76 TeV are 83% pions, 13% kaons, and 
4% protons, with the pion fraction increasing to 84% for 
more peripheral (50-60%) collisions. Hence, I use an ac¬ 
cordingly averaged effective mass for m as described in 
detail in Q- 

Since only the ratio mjpT enters the Jacobian, one 
can also fix the mass at the pion mass m = and cal¬ 
culate the corresponding effective transverse momentum 

from < pT,eff >= '^ 7 TJy=o/'^'^ — Jy=o with the exper¬ 
imentally determined Jacobian Jy=o at 90^, see 0 for 
2.76 TeV PbPb. The values of mjpT used in the calcu¬ 
lations shown in Fig. [5] are mjpT = 0.466; 0.349; 0.585 
for ^/snn = 0.13; 0 . 2 ; 2.76 TeV, respectively. 

The above expression for the transformation from y to 
77 has the limits 7 / ^ 77 — ln(m/pT) for m « pt, and 
7 / ^ 77 for pt « m. Since most of the produced charged 
hadrons at a LHG energy of 2.76 TeV are pions, the limit 
77 ~ 77 at small transverse momenta - very forward an¬ 
gles - is reached for charged hadrons at larger values of 
77 than for protons (net protons determine the value of 
the beam rapidity). Hence, the dN/dy distribution for 
charged hadrons which are mostly pions can easily extend 
beyond 7 /b earn* 


V. CONCLUSION 

The analysis of the energy dependence of charged- 
hadron pseudorapidity distributions in AuAu collisions 
at RHIC energies, and PbPb collisions at LHC ener¬ 
gies in the phenomenological nonequilibrium-statistical 
relativistic diffusion model reveals the expected ln(s)- 
dependence for the total particle content of the two frag¬ 
mentation sources, but a In^s-dependence for the total 
charged hadron content of the gluon-gluon source. Mod¬ 
ifying the conclusion of an initial investigation 0 , it is 
only in a limited energy region of about 100 GeV to 2.76 
TeV that this dependence may be approximated by a 
power law. 
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